arXiv:1505.02165vl [astro-ph.CO] 8 May 2015 


Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 


Printed 12 May 2015 (MN MbX style file v2.2) 


Next Generation Cosmology: Constraints from the Euclid Galaxy 
Cluster Survey 

B. Sartoris^’^, A. Biviano^, C. Fedeli^’"^, J. G. Bartlett^, S. Borgani^’^’^, M. Costanzi^, 

C. L. Moscardini^’^’'^, J. Weller^’^®’^\ B. Ascaso^^, S. Bardelli^, 

S. Maurogordato^^, and P. T. P. Viana^"^’^^ 

* Dipartimento di Fisica, Sezione di Astronomia, Universita di Trieste, Via Tiepolo 11,1-34143 Trieste, Italy 
^ INAF/Osservatorio Astronomico di Trieste, Via Tiepolo 11,1-34143 Trieste, Italy 
^ INAF/Osservatorio Astronomico di Bologna, Via Ranzani 1,1-40127 Bologna, Italy 
^ INFN, Sezione di Bologna, Viale Berti Pichat 612,1-40127 Bologna, Italy 
^ APC, AstroParticule et Cosmologie, Universite Paris Diderot, CNRS/1N2P3, CEA/lrfu, 

Observatoire de Paris, Sorbonne Paris Cite, 10, rue Alice Domon et Leonie Duquet, 75205 Paris Cedex 13, France 
^ INFN, Sezione di Trieste, Via Valerio 2, 1-34127 Trieste, Italy 

^ Universitdts-Sternwarte Munchen, Fakultdt fur Physik, Ludwig-Maximilians Universitdt Munchen, Scheinerstr. 1, D-81679 Muiichen, Germany 
® Dipartimento di Fisica e Astronomia, Alma Mater Studiorum Universita di Bologna, viale Berti Pichat, 612, 40127 Bologna, Italy 
^ Aix Marseille Universite, CNRS, LAM (Laboratoire dAstrophysique de Marseille) UMR 7326, 13388, Marseille, France 
Excellence Cluster Universe, Boltzmannstr. 2, D-85748 Garching, Germany 
** Max Planck Institute for Extraterrestrial Physics, Giessenbachstr. 1, D-85748 Garching, Germany 

GEPl, Observatoire de Paris, CNRS, Universite Paris Diderot, 61, Avenue de VObservatoire 75014, Paris Prance 

Laboratoire J.-L. LAGRANGE, CNRS/UMR 7293, Observatoire de la Cote dAzur, Universite de Nice Sophia-Antipolis, 06304 Nice Cedex 4, Prance 

Instituto de Astrofi'sica e Ciencias do Espaqo, Universidade do Porto, CAUP, Rua das Estrelas, 4150-762 Porto, Portugal 

Departamento de Fisica e Astronomia, Faculdade de Ciencias, Universidade do Porto, Rua do Campo Alegre 687, 4169-007 Porto, Portugal 


12 May 2015 


ABSTRACT 

We study the characteristics of the galaxy cluster samples expected from the European Space 
Agency’s Euclid satellite and forecast constraints on parameters describing a variety of cos¬ 
mological models. The method used in this paper, based on the Fisher Matri x approach, is 
the sa me one used to provide the constraints presented in the Euclid Red Book (iLaureiis et alJ 
1201 ih . We describe the analytical approach to compute the selection function of the photo¬ 
metric and spectroscopic cluster surveys. Based on the photometric selection function, we 
forecast the constraints on a number of cosmological parameter sets corresponding to dif¬ 
ferent extensions of the standard ACDM model, including a redshift-dependent Equation of 
State for Dark Energy, primordial non-Gaussianity, modified gravity and non-vanishing neu¬ 
trino masses. Our results show that Euclid clusters will be extremely powerful in constraining 
the amplitude of the matter power spectrum erg and the mass density parameter Qm. The dy¬ 
namical evolution of dark energy will be constrained to Awo = 0.03 and Awa = 0.2 with 
free curvature Qj:, resulting in a (wo, Wa) Figure of Merit (FoM) of 291. Including the Planck 
CMB covariance matrix, thereby information on the geometry of the universe, improves the 
constraints to Awo = 0.02, Avv^ = 0.07 and a FoM= 802. The amplitude of primordial non- 
Gaussianity, parametrised by /nl, will be constrained to A/nl - 6.6 for the local shape sce¬ 
nario, from Euclid clusters alone. Using only Euclid clusters, the growth factor parameter 
y, which signals deviations from General Relativity, will be constrained to Ay = 0.02, and 
the neutrino density parameter to AQy = 0.0013 (or A2»*y = 0.01). We emphasise that 
knowledge of the observable-mass scaling relation will be crucial to constrain cosmological 
parameters from a cluster catalogue. The Euclid mission will have a clear advantage in this 
respect, thanks to its imaging and spectroscopic capabilities that will enable internal mass 
calibration from weak lensing and the dynamics of cluster galaxies. This information will 
be further complemented by wide-area multi-wavelength external cluster surveys that will 
already be available when Euclid flies. 
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1 INTRODUCTION 

According to the hierarchical scenario for the formation of cos¬ 
mic structures, galaxy clusters are the latest objects to have formed 
from the collapse of high density flu ctuations filtered on a typi - 
cal scale of ~ 10 comoving Mpc (e.e. lKravtsov & Borganill2012h . 
Since galaxy clusters provide information on the growth history 
of structures an d on the underlying cosmolog ical model in many 
ways (see, e.e.. lAllen. Evrard & Manta 1201 ih . they have played 
an important role in delineating the current standard ACDM cos¬ 
mological model. As a matter of fact, the number counts and 
spatial distribution of these objects have a strong dependence on 
a number of cosmological parameters, especially the amplitude 
of the mass power spectrum and the matter content of the Uni¬ 
verse. The evolution with redshift of the cluster number den¬ 
sity and correlation function can be employed to break the de¬ 
generacy between these two parameters, and thus can provide 
constraints on the cold Dark Matter ( DM henceforth) and Dark 
Energy (DE) density parameter s (e.g.. IWang & Steinhardtl 


Haimaiij^Mohr&Jfolde^ 2001]^ Wellerj_Battvej&Knei^ JzytM 

Battve & Welle J2OO3I: I Allen. Evrard & Mantj|201ll:ISartoris et alj 


199S 


2002; 


2012h. Eurthermore. a numb er of studies (e.g.. lCarbone et al .I2OI2I; 


Costanzi et alj|20133 . l2014h have also shown that clusters can be 


used to constrain neutrino properties, because massive neutrinos 
would directly influence the growth of cosmic structure, by sup¬ 
pressing the matter power spectrum on small scales. More gener¬ 
ally, since the evolution of the cluster population traces the growth 
rate of density perturbations, large surveys of clusters extending 
over a wide redshift interval have the potential of providing strin¬ 
gent constraints on any cosmological model whose deviation from 
ACDM leaves its imprint on this growth. 

Over the past decade, surveys of galaxy clusters for cos¬ 
mological use have been constructed and analys ed, based on 
observations at different wavelengths: X-ray (e.g. iBorgani et alj 
l200ll;IVikhlinin et al .11200^; IClerc et alJl2012tlRaDetti et al.ll2013h 


sub-mm, throug h the ISun yaev & Zeldovichl jl972|) dis tortion 
(SZ henceforth, |Stanisze wski_eL^ |2009|; LBenspn_eta D I2013I: 


IPlanck Collaboration et al.l l2014bl ; iBurenin & VikhlininI EoTll . 

and optical jRozo et al.l 1201 ol) bands. Eurther improvements 
can be obtained from the s patial cluster i ng of galaxy clus - 
ters dSchuecker et a 0 I2OO3I; iHiitsil l2010l: iMana et alj l2013h . 
The resulting cosmological constraints turn out to be com¬ 
plementary to those of ot her cosmological probes such as 
type la supemovae (e.g. 


. her cosmotogicat p r 

iBetoule et al. I20141). 


Cosmic Mi- 


Background (CMB) radi ation te.g.. lHinshaw et alJl201.3l: 


crowave _ 

IPlanck Collaborat ion et al.ll2014ah. the B aryon Acoustic Oscilla- 
tions (BAOs; e.g., Anderson et al.|[2014) . and cosmic shear (e.g. 
iKitching et alJ2of4 ). These cluster catalogues are however charac¬ 
terised either by a large number of objects that cover a relatively 
small redshift range, or rather small samples that span a wide red¬ 
shift range. Ideally, in order to exploit the redshift leverage with 
good statistics, one should have access to a survey that can provide 
a high number of well characterised clusters over a wide redshift 
range. 

One future mission that will achieve this goal will be the 
European Space Ag ency (ESA) Cosmic Vision mission Eucli^ 
jUaureiis et alj|20f]|) . Planned for launch in the year 2020, Euclid 
will study the evolution of the cosmic web up to redshift z ~ 2. 
Although the experiment is optimised for the measurement of cos¬ 
mological Weak Lensing (WL, or cosmic shear) and the galaxy 

* http://www.euclid-ec.org 


clustering, Euclid will also provide data usable for other important 
complementary cosmological probes, such as galaxy clusters. Clus¬ 
ter detection will be possible in three different ways: i) using pho¬ 
tometric data; ii) using spectroscopic data; and Hi) through gravi¬ 
tational (mostly weak) lensing, which may be combined for more 
efficiency. In this paper, we will perform our analyses by using the 
photometric cluster survey (see Section|3, where the cluster detec¬ 
tion method is not dissimilar from t hat used to detect low-redshift 
SDSS clusters dKoester et alj|20(y^ . However, thanks to the use 
of Near Infrared (NIR) bands, Euclid will be capable of detecting 
clusters at much higher redshifts (z ~ 2) over a similarly large area. 
The sky coverage of Euclid will reach 15,000deg^, almost the en¬ 
tire extragalactic celestial sphere. The characteristics of the Euclid 
spectroscopic survey and its possible use for the calibration of the 
mass-observable relation will be discussed in Appendices A and B, 
respectively. 

One fundamental step for the cosmological exploitation of 
galaxy clusters is the definition of the relation between the 
mass of the host DM halo and a suitable observ able quan¬ 
tity (e.g.. lAndreon & Hurnll2012l ; lOiodini et al.ll201^ . Many ef¬ 
forts have been devoted to the calibration of the observable- 


20 id; Planck Collaboration et al.1 

2011; Reichert et al. 2011; 

Rozo et al.ll201ll; Rvkoff et alj2012l; 

Ettoril2013l;lRozo et alj2014 ; 


are at the same time precise (i.e. characterised by a small scat¬ 
ter in the scaling against cluster mass) and robust (i.e. rel- 
atively insensitive to the details of cluster astrophysics) (e.g. 
iKravtsov. Vikhlinin & Nag^l2006h . In the case of Euclid, an in¬ 
ternal mass calibration will be performed through the exploitation 
of spectroscopic and WL data of the wide Euclid survey (see Ap¬ 
pendix!^, of tho deep Euclid survey of 40deg^, 2 magnitudes 
deeper than the wide survey. The deep survey will be particularly 
useful in adding constraints on the evolution of the observable-mass 
scaling relation at z > 1. 

These Euclid internal data will provide a precise calibra¬ 
tion of the relation between cluster richness, which characterises 
photometrically-identified clusters, and their actual mass. Eur¬ 
thermore, it will be possible to cross-correlate Eu clid data with 
data from o ther cluster surveys - such as eRosita dMerloni et alj 
2012h. XCS (iMehrtens et alj2012l) . the South Pole Telescope (SPT, 


Carlstrom et al 


HF 

(ACT. lMarriage et alj|20l7h - to further improve the mass calibra¬ 
tion of Euclid clusters. 

The aim of this paper is to forecast the strength and the 
peculiarity of the Euclid cluster sample in constraining the pa¬ 
rameters describing different classes of cosmological models that 
deviate from the concordance ACDM paradigm. We first con¬ 
sider the case of a dynamical evolution of the DE compo- 
ne nt, using the two-parameter fu nction al form origina lly proposed 
by IChevallier & Polarskil ( 1200 ih and iLindeJ (l2()03h . The same 
parametrisati on has been used in the Da rk Energy Task Eorce re¬ 
ports IDETE: [Albrecht et alJl2006Ll2009h to estimate the constrain¬ 
ing power of different cosmological experiments. Second, we al¬ 
low for the primordial mass density perturbations to have a non- 
Gaussian distribution. Third, we explore the effect of deviations 
from General Relativity (GR) on the linear growth of density per¬ 
turbations. Einally, we consider the case of including massive stan¬ 
dard neutrinos. 

The structure of this paper is the following. In Section [2] we 
describe the approach used to estimate the Euclid cluster selection 
function of the photometric survey. In Section we describe the 


201 Ih. and the Atacama Cosmology Telescope 


© 0000 RAS, MNRAS 000, 000-000 


























































































































Next Generation Cosmology: Constraints from the Euclid Galaxy Cluster Survey 3 



redshift 

Figure 1. Number Nsoo,c of cluster galaxies within rsoo,c (black curves), and 
ScTgeid where CTfigu is the rms of the field counts within the same radius, and 
within the adopted 3Azp cut (red curves). These counts are shown down to 
the limiting magnitude of the Euclid survey, Hab = 24, as a function of 
redshift for clusters of different masses, log(M 2 oo,c/MQ) = 14.5,14.0,13.5 
(solid, dot-dashed, dashed line, respectively), where masses are defined with 
a mean overdensity of 200 times the critical density of the universe at the 
cluster redshifts. 

Fisher Matrix approach used to derive constraints from the Euclid 
cluster survey on cosmological parameters. In Section|4l we briefly 
describe the characteristics of the different cosmological models 
we consider. In Section [S] we show our results on the number of 
clusters that the wide Euclid survey is expected to detect as a func¬ 
tion of redshift and the constraints that will be obtained on the cos¬ 
mological parameters using the cluster number density and power 
spectrum. Finally, we provide our discussion and conclusions in 
Section[6l We present the analytical derivation of the spectroscopic 
selection function in Appendix IaI and the calibration of the cluster 
observable-mass relation in Appendix IB] 


2 GALAXY CLUSTER SELECTION IN THE Euclid 
PHOTOMETRIC SURVEY 

In this Section, we adopt the cosmol ogical parameter values of 
the con cordance ACDM model from IPlanck C^laboration et alj 
( l2014ah . Ho = 67 km s^'Mpc^' for the Hubble constant, On, = 0.32 
for the present-day matter density parameter, and Oj. = 0 for the 
curvature parameter. 

To determine the selection function of galaxy clusters in 
the Euclid photometric survey, we adopt a phenomenological ap¬ 
proach. We start by adopting an average universal luminosity 
functio n (LF hereafter) for cluster galaxies. lUin, Mohr & Stanford 
( l2003h evaluated the Xj-band LFs of cluster galaxies out to a radius 
fsoo.c for several nearby clusters. The radius is defined as the ra¬ 
dius of the sphere that encloses an average mass density A times the 
critical density of the Universe at the cluster redshift. These clus¬ 
ter L Fs were parametrised using Schechter functions l lSchechteil 
Il976h . We adopt the averages of the nor malisations and char¬ 
acteris tic luminosities listed in Table 1 of lUin. Mohr & Stanford 
( l2003h for the 27 nearby clusters included in that analysis, corre- 
sponding to <^* = 6.4 Mpc~^ and M* = -24.85. Also, following 
iLin. Mohr & Stanford l l2003h . we use a faint-end slope a = -1.1, 
as confirmed in the r-band deep spect roscopic analysis of two 
nearby clusters bv iRines & Gelleil bOOSh . 


Concerning the behaviour of the cluster LF at z > 0, 
there is no conclusive observational evid ence on the evolution 
of the LF faint-end s l ope p arameter a dMancone et al.l 12012l : 
IStefanon & Marchesinil l2013h . Therefore, we assume it to be 
redshift-invariant. Also, the observed constancy of the richness 
vs. mass relation for clusters up to z ^ 0.9 dLin et al.l l2006l : 
IPoggianti et all l20ld : lAndreon & Congd^ Eoidh suggests that 
there is no redshift evolution of (f>*, apart from the cosmological 
evolution of the critical density, which scales as H‘(z). 

We assume the M* parameter to change with z ac- 
cording to passive evolu tion models of stellar populations 
dKodama & Arimotall997l) . This assumption is justified because 
emission in the X, band is not strongly influenced by young stel¬ 
lar ge nerations, and it is supported by observations dMancone et alj 
I2OI2L and references therein), at least for clusters more massive 
than ~ lO*'^ Mq. For clusters of lower mass, some high-z surveys 
have found evidence for deviation from passive evolution of M* 
dMancone et al .l201(]l:lTran et alj20ld : lBrodwin et al.l2^ . How¬ 
ever, the current observational evidence does not allow us to pre¬ 
cisely parametrise M* evolution to z > 1 and low cluster masses, 
and we prefer to keep our conservative assumption of passive evo¬ 
lution over the full cluster mass range. _ 

W e apply the early-type ^-correction of iMannucci et alj 
d200lh to the M* magnitudes. This correction should be the most 
appropriate for galaxies in clust ers, which are mostly early-type 
even at relatively high redshifts dPostman et alj|200l : ISmith et alj 
l2005h . We finally convert the Ks magnitudes into the Euclid band 
Hab using the mean rest-frame colour for cl uster galaxies, H - 
Ks = 0.26 (we average the values provide d bv iBoselli et alj[l997l : 
Ide Pronris et al.ll998l : lRamella et al.l2004h . an d adopting the trans¬ 
formation to the AB-system Hab = H + 1.37 dCiliegi et al.ll2005]) . 
We thus obtain the cluster LFs in the Hab band at different redshifts. 

By integrating these LFs down to the apparent magnitude 
limit of the wide Euclid photometric survey (Hab — 24, see 
iLaureiis et al.ll201ll) . we then evaluate nsoo.c, namely the redshift- 
dependent number density of cluster galaxies within rsoo.c- The 
number of cluster galaxies contained within a sphere of radius 
rsoo.c (i-e. the cluster richness) is then Ni(ya,c = = 

8/3;rn5oo,cGM5oo,c/[500H^(z)], where the last equivalence follows 
from the relation between rjoo.c and Mjoo.c, the mass within a mean 
overdensity of 500 times the critical density of the universe at the 
cluster redshifts. Note that the dependence of Al5oo,c on H^^(z) is 
only apparent, since (p*, and hence nsw.c, scales as H^(z)- The z- 
dependence only comes in as a result of the fixed magnitude limit 
of the survey and the passive evolution of galaxies. In Fig. [T] we 
show Al5oo,c(z) for clusters of three different masses: logM2oo,c = 
13.5, 14.0, and 14.5 (black curves). To convert from t o M 200 .C 
we adopt a NFW profile jNavarro. Frenk & White|[l9^ with a 
m ass- and redshift-dep endent concentration given by the relation 
of iDe Boni et alj ( 12013l . 2"“^ relation from top in their Table 5). 

We then estimate the contamination by field galaxies in 
the cluster area. We take the estimate of the number density of 
field galaxies down to Hab = 24 from the H-band counts of 
iMetcalfe et all ( I 2 OO 6 L see their Table 3), nfldd - 33 arcmin 
an esti mate that is in agree ment with the Euclid survey require¬ 
ments dLaureiis et al.ll201]I) . Multiplying this density by the area 
subtended by a galaxy cluster at any given redshift we obtain the 
number of field galaxies that contaminate the cluster field-of-view, 
Afield = «fieid ?tr500,c’ where rjoo.e is in arcmin. 

The number of field contaminants can be greatly reduced by 
using photometric redshifts, Zp. These will be obtained to the re¬ 
quired accuracy of Azp = 0.05(1-1-Zc), by combining the Euc/idpho- 
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redshift 


M2oo,c ~ 8 X 10*^ Mq roughly corresponds to ~ 80% completeness 
at all redshifts z < 2. 

The shape of the selection functions shown in Fig.[2]is some¬ 
what counter-intuitive because it is higher at z ~ 0.2 than at z ~ 0.7. 
Naively one would expect that clusters of lower mass would be eas¬ 
ier to detect at lower redshifts. We find that this shape is related to 
the relative importance of cosmic variance and Poisson noise in the 
contaminating field counts. Cosmic variance drives the shape of the 
selection function at z < 0.5 and Poisson noise at higher redshifts. 
If we select clusters at a higher overdensity (e.g. A^. = 2500 rather 
than Ac = 500), the relative importance of cosmic variance and 
Poisson noise changes in a way to flatten the selection function at 
z < 0.5. In reality, observers do not select clusters at given A^, so 
our estimate of the selection function must be considered only as an 
approximation. At the end of Section we comment on the effect 
of taking a flat selection function out to z = 2. 


Figure 2. Galaxy cluster mass selection function for the Euclid photometric 
survey. Solid and dashed lines are for detection thresholds Alsoo.c/cfieid = 3 
and 5, respectively. 


tomet ric survey with auxiliary ground-based data dLaureiis et alJ 
l201lh . One can safely consider as non-cluster members all those 
galaxies that are more than 3Azp away from the mean cluster red- 
shift Zc. The mean cluster redshift will be evaluated by averaging 
the photometric redshifts of galaxies in the cluster region, and addi¬ 
tionally including the (few) spectroscopic galaxy redshifts provided 
by the Euclid spectroscopic survey (see Appendix lAt. 

In order to determine the fraction of field galaxies, /(Zc), with 
photometric redshift Zp in the range ±3 x 0.05(1 -f Zc) at any given 
Zc, we need to estimate the photometric redshift distribution of an 

= 24 limited field survey. To this aim we consider the photo¬ 
metric re dshift distribution of galaxies with Hab ^ 24 in the cat¬ 
alogue of lYang ai] ( l2014h . We find /(zc) = 0.07,0.23,0.34, and 
0.33 at Zc = 0.2,0.8,1.4, and 2.0, respectively. 

Finally, we evaluate the rms, crgeid, of the field galaxy counts 
f{Zc)Nfic\i^ by taking into account both Poisson noise and cos¬ 
mic variance. For the latter we use the IDL code quickcv of 
John Moustakaj^ for cosmic variance calculation. In Fig. [T] we 
show 3crfleid as a function of redshift, in clusters of log M 200 .C - 


13.5,14.0, and 14.5. 

The ratio between the cluster galaxy number counts and the 
field nns, Nsoq.c/o'rm^ gives the significance of the detection for a 
given cluster. The cluster selection function is the limiting cluster 
mass as a function of redshift for a given detection threshold. This is 
shown in Fig.[2]for two thresholds, Ajoo.c/crflcid = 3, and 5. This se¬ 
lection function is only mildly dependent on redshift. The limiting 
cluster mass for the lowest selection threshold (Agoo.c/o'fieid = 3) 
is M2oo,c ~ 8 X 10*^ Mq, lower than the typical mass of richness 
class 0 clusters in th e lAbell. Corwin & OlowinI ( Il989h catalogue 
dPopesso et al.ll2012l) . It is also similar to the limiting mass of the 
selection function of SPSS clusters ide ntified by the maxBCG al¬ 
gorithm (see Fig. 3 inlRozo_et_alj201^, and to the typical mass of 
the clusters identified by Brodwin et alj (l2007h up to z ~ 1.5 using 
Zp in an IR-selected galaxy catalogue. Preliminary tests based on 
running cluster finders on Euclid mock^ show that the mass limit 


^ https://code.google.eom/p/idl-moustakas/source/browse/trunk/impro/cosmo/ 
quickcv.pro?r=617 

^ http://wiki.cosmos.esa.int/euclid/index.php/EC_SGS_OUXE3. Access 
restricted to members of the Euclid Consortium. 


3 FISHER MATRIX ANALYSIS 


Before presenting our forecasts for the cosmological constraints we 
now briefly describe the Fisher Matrix (FM hereafter) formalism 
that we use to derive these constraints. 

The FM formalism is a Gaussian approximation of the likeli¬ 
hood around the maximum to second order and it is an efficient way 
to study the accuracy of the estimation of a vector of parameters p 
by using independent data sets. The FM is defined as 


Fa0 = - 


l d^\n£ \ 

\ dpadp/i I ’ 


( 1 ) 


where £ is the likelihood of the observable quantity te.g. lPodelsonl 
l2003h . In our FM analysis we combine three different pieces 
of information: the galaxy cluster number density, the clus¬ 
ter power spectrum, and the prior knowledge of cosmologi- 
cal parameters as derived from the Planck CMB experiment 
jPlanck Collaboration et al.ll2014al) . To quantify the constraining 
power of a given cosmological probe on a pair of j oint parameters 
(p,-. P i) we use the Figure of Merit (FoM henceforth: I Albrecht et alj 
1200 ^ 


FoM = - 

,^det [Cov(p„Pj)] 


( 2 ) 


where Cov(pi,pj) is the covariance matrix between the two param¬ 
eters. With this definition, the FoM is proportional to the inverse 
of the area encompassed by the ellipse representing the 68 per cent 
confidence level (c.l.) for mod el exclusion. 

As des cribed in detail in|Sartgris_gtal 


_ (l2010h . we follow the 

approach of lHolder. Haiman & Mohi ( 2001 ) and define the FM for 


the cluster number counts as 


_ ^ y dNip,dNfp, 1 


In the previous equation, the sums over i and m run over redshift 
and mass intervals, respectively. The quantity is the number of 
clusters expected in a survey with a sky coverage Psky, within the 
Ath redshift bin and m-th bin i n observed mass M°'°. This can be 
calculated as dLima & Htill2005h 

Psky 


N, 


£,m 


, dV 


X 

dM°'^ 

- dMn{M,z)p{M° 

Jm* Jo 


W), 


(4) 
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where dVIdz is the cosmology-dependent comoving volume ele¬ 
ment per unit redshift interval. The lower observed mass bin is 
bound by = Mthr(z), where Mthiiz) is defined as the threshold 

value of the observed mass for a cluster to be included in the sur¬ 
vey (see Fig. 13. For the halo mass func tion n{M, z) in equat ion (|4l(, 
we assume the expression provided by [Tinker et all ( l2008h . Since 
the Euclid selection function has been computed for masses at 
= 200 with respect to the critical density, we use the [Tinker et all 
( [2008h mass function parameters relevant for an overdensity of 
Am- = 200/f2„,(z) with respect to the background density. We note 
that in equation @ we have implicitly assumed that the survey sky 
coverage is independent of the observed mass, which may not 
necessarily be the case if the sensitivity is not constant over the 
survey area. 

In equation (|4j, is the probability to assign an ob- 

served mass M°'° t o a galaxy cluster with true mass M. Following 
[Lima & Hij ( [2005h . we use a lognormal probability density, namely 


where 


p(M°''\M) = 


exp[-x2(M“'’)] 


^ lnM°^-lnMh.a.-lnM 


(5) 


( 6 ) 


In the above equation InMbias is the bias in the mass estimation, 
which encodes any scaling relation between observable and true 
mass and should not be confused with the bias in the galaxy distri¬ 
bution. CTinM is the intrinsic scatter in the relation between true and 
observed mass (see Section^. By inserting equation Q into equa¬ 
tion we obtain the expression for the cluster number counts 
within a given mass and redshift bin. 


Nt.^ 


flaky C'*' , dV 

8tr X ^ * 



dMn(M,z) [erfc(x„,) - erfc(ji;,„+i)] , 


(7) 


where erfc(x) is the complementary error function and Xm = 



The FM for the averaged redshift-space cluster power spec¬ 
trum within the f’-th redshift bin, the /n-th wavenumber bin, and the 
i-th angular bin can be written as 


„ 1 ^ d\nP(pi,lc„,zt) d\nP(jj.i,k„„zt) ,y.,. 


<9p« 


(8) 

(e.g.. [Tegmark[[l99'^ : [Feldman. Kaiser & Peacock! 1994) . where the 
sums over £, m, i run over bins in redshift, wavenumber, and cosine 
of the angle between k and the line of sight direction, respectively. 
The quantity V‘^*^(p,, k,„,zi) represents the effective volume accessi¬ 
ble to the survey at redshift zt and wavenumber k iTegmark[[l997l : 
[Sartoris et alj[201^ . and reads 


= V'o(Zf) 


h{ze)P(piNm, Ze) 

1 + h(zt)P{piNm,zt) ' 


(9) 


In the above equation, Vf,{zi) is the total comoving volume con¬ 
tained in the unity redshift interval around zt, while h{z() is the av¬ 
erage number density of objects included in the survey at redshift 
Zt, 


nizt) = 



dM n{M, Zt) erfc {x[Mthr(zr)l . 


( 10 ) 


The cluster power spectrum averaged over the redshift bin, appear¬ 
ing in equation can be written as 

1 n*' dV , - 

F(p„ k„„ Zt) = — dz — n^iz) Pijii, km, z) , (11) 

S t Jzi dz 

where the normalisation factor S i reads 

dV , 

St= dz — h\z). (12) 

Jzt dz 

[Sartoris et alj ( [2012h pointed out the importance of taking into 
account the contribution of cluster redshift space distort i ons fo r 
constraining cosmological parameters. Following [Raise I ( [ 19871) . 
we calculate the redshift-space cluster power spectrum P(pi, k^, Zt) 
in the linear regime according to 

Hpi,k,„,zt) = [heizt) + f(zt)N^ Fl (F„„ Zt) , (13) 

where the power spectrum acquires a dependence on the cosine p of 
the angle between the wavevector k and the line-of-sight direction. 
In the above equation, b^gizi) is the linear bias wei ghted by the 
mass function (see equation 20 in [Sartoris et al.[|^10l) . 


bcsizt) = 


— r 

n{zt) Jo 


dMn(M,Zf)erfc{x[M,hr(Z(’)]| b(M,zt). (14) 


The function f(a) = dlnD(a)/dlna is the logarithmic derivative 
of the linear growth rate of density perturbations, D(a), with re¬ 
spect to the expansion factor a. Piik,„,zt) is the linear matter power 
spectrum in real space, that we c alculate using the CLASS code 
( [Bias. Lesgourgues & Tramir201lh . For the DM halo bias b{M,z) 
we use the expression provided by [Tinker et alj ( [201(l[) . 

Both the power spectrum and the number counts FMs (equa¬ 
tions [3 and [3 are computed in the redshift range defined by the 
Euclid photometric selection function shown in Fig. [3 namely 
0.2 < z < 2, with redshift bins of constant width Az = 0.1. 
We note that the limiting precision with which the redshift Zc 
of a cluster is determined in the photometric survey is given by 
0.05(1 -I- Zc)/A^5gQ^, where Ajoo.c is the total number of galaxies as¬ 
signed to the cluster. Therefore, the bin width is always larger than 
the largest error on redshift expected from the Euclid photometric 
survey (see Section|3- In equation l[3, the observed mass range ex¬ 
tends from the lowest mass limit determined by the photometric se¬ 
lection function (Tfthr(z), see Fig.|3 up to loglMob/Mg) < 16, with 
Alog(Mob/Afo) = 0.2. In the computation of the power spectrum 
FM (equationlDl, we adopt k^ax = 0.14Mpc“*, with Alog(kMpc) = 
0.1. Finally, the cosine of the angle between k and the line of sight 
direction , p, runs in the rang e -1 ^ p < 1 with 9 equally spaced 
bins (see Sartoris et alj[201 


4 COSMOLOGICAL AND NUISANCE PARAMETERS 

In this Section we discuss the cosmological parameters that have 
been included in the FM analysis in order to predict the constrain¬ 
ing power of the Euclid photometric cluster survey and we describe 
the peculiarity of all the analysed models. As a starting point, we 
consider all the standard cosmological parameters for the concor- 
da nce ACDM model, whose fiducial v alues are chosen by follow¬ 
ing [pianck]Con^gr^ion|er^ | [2014ah : Dm = 0.32 for the present- 
day total matter density parameter, erg = 0.83 for the normalisation 
of the linear power spectrum of density perturbations, Db = 0.049 
for the baryon density parameter, Hq = 67kms“'Mpc“’ for the 
Hubble constant, and ng = 0.96 for the primordial scalar spectral 
index. We also allow for a variation of the curvature parameter, 
whose fiducial value Oj. = 0 corresponds to spatial flatness. 
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4.1 Model with dynamical Dark Energy 

In addition to the ACDM parameters, we also include parame¬ 
ters describing a dynamical evolution of the DE component. In 
the literature there are a number of models, characterised by dif¬ 
ferent parametrisation of the DE Equa tion of State (EoS hence¬ 
forth) evolution (e.g., IWettericlJ |2004) . In this paper we study 
the pa ram etrisation origi nally proposed by Ichevallier & Polarskil 
(l200lh and lLind^ l l2003h and then adopted in the DETE. We label 
this parametrisation as the CPL DE model, according to which the 
DE EoS can be written as 


w(a) = Wo + Wa(l - a) . (15) 

We use Wo = -1 and Wa = 0 as reference values for the two model 
parameters. Thus, the cosmological parameter vector that we use in 
this first part of our EM analysis reads 


p = {Dn„a-8, Wo, Wa,Dt,Db,^fo,nsl • (16) 


The constraints on the DE dynamical evolution obtained by com¬ 
bining Planck C MB data with WMAP po l arisatio n and with 
ESS information jPlanck Collaboration et al.l 12014ah . are wo = 
-1.04);ggg and Wa S 1.3 (95 per cent c.l.) assuming Slu = 
0. Currently, the evolution of the cluster number counts alone 
does not constrain the DE equation of state parameters. However, 
iMantz et al.l ( 12014h were able to obtain: wo = -1.03 ± 0.18 and 
Wa = -0. Dq® (assuming O*. = 0), by using CMB power spec¬ 
tra (1-year Planck data, SPT, ACT), SNIa, and BAO data at differ- 
ent red shifts (plus WMAP polarisation: IPlanck Collaboration et alj 
l2014ah . 


Despite these weak c onstraints on the CPL DE parametrisa¬ 
tion dVikhlinin et alj|2009l) . cluster counts are powerful probes of 
the amplitude of the matter power spectrum. For instance, erg is 
constrained at th e level of ~ 8 per c ent both with optically selected 
SDSS clusters ( Rozoet_^ 120101) . and with SZ selected SPT 
clusters dBenson et al.ll201^ . Moreover, clusters help breaking the 
degeneracy between erg and Dm in CMB datasets, improving the 
constraints on the amplitude of the matter power spectrum by a fac - 
tor of ~ 2 with respect to CMB constraints alone dRozo et al.l2010l) . 


4.2 Model with primordial Non-Gaussianity 

We extend the standard cosmological model by allowing primor¬ 
dial density fluctuations to follow a non-Gaussian distributio n 
(e.e.. lBartolo et al]|2004l : lDesiacaues & Selia^l2010l : IWang|l20l^ . 

When this happens, the distribution of primordial fluctuations in 
Bardeen’s gauge-invariant potential O cannot be fully described 
by a power spectrum - commonly parametrised by a power-law, 
P<t{k) = (where k = ||k||) - rather we need higher-order 

statistics such as the bispectrum 6o(ki,^ 2 ,Different models 
of inflation are known to produce different shapes of this bispec¬ 
trum. Here we consider only the so-called local shape, where the 
bispectrum strength is maximised for squeezed configurations, in 
which one of the three momenta kj is much smaller than the other 
two. 

Within the local shape scenario, we adopt the commonly used 
way to parametrise the primordial non-Gaussianity, which allows 
us to write Bardeen’s gauge invariant potential as the sum of a lin¬ 
ear Gaussian term and a non-linear second- order term that encapsu¬ 
lates the deviation from G aussianity (e.e.. lSaloDek & Bondlll990l : 
iKomatsu & Spergell200lh : 

O = Og -t /nl (®| - <0|>) , (17) 


where the free dimensionless parameter /nl parametrises the 
deviation from the standard Gaussian scenario. We stress that 
there is some ambiguity in the normalisation of equation 03- 
We adopt t he ESS convention (as opposed to the CMB con¬ 
vention, see Pillepich. Porciani & HahnlBoid : iGrossi et^ l2007l : 


ICarbone. Verde & Matarresell2008ah where O is linearly extrapo¬ 
lated at z = 0 for defining the parameter /nl- The relation between 
the two normalisations is /nl = D{z = oo)(l -i- z)/n]^®/D(z = 0) - 
1.3/™®, where D{z) is the linear growth factor with respect to the 
Einstein-de Sitter cosmology. 

If the density perturbation field is non-Gaussian and has a pos¬ 
itively (negatively) skewed distribution, the probability of forming 
large overdensities - and thus large collapsed structures - is en¬ 
hanced (suppressed). Thus, the shape and the evolution of the mass 
funct i on of DM halos change (e.g., Matairese. Verde & Jimene3 
l200(llGrossi et alj|2009l: ILo V erde et alJl26oik Following the pre¬ 
scription by LoVerde et alj ( l2008h one can modify the mass func¬ 
tion n{M, z) in equation 0 to take into account the non-Gaussian 
correction as follows 


n(M,z) = n^^\M,z) 


^ nps{M,z) 
nfJ(M,z) 


(18) 


In the previous equation, z) is the mass function in the refer- 

ence Gaussian model, wh ile nps{M,z) and «pg'(M,z) represent the 
IPress & Schechteil (Il974h mass functions in the non-Gaussian and 
reference Gaussian models, respectively (see the full equations in 
ISartoris et alj|201^ . 

In non-Gaussian scenarios the large-scale clustering of DM 
halos also changes. This modification is quite important because it 
alters in a fairly unique way the spatial distrib ution of tracers of 
the cosmic structure, including galaxy clusters lDalalrt ^l2008l : 


iMatarrese & Verd^l2008l : iGiannantonio & Porcianill201(]ll . Specifi 

cally, the linear bias acquires an extra scale de pendence due to pri¬ 
mord ial non-Gaussianity, and can be written as dMatarrese & Verd3 

I 2 OO 8 I 1 

KM, z, k) = K°\M, z) -t [K°\M, z) - 1 ] <5e(z)r«(k) , (19) 

where TKk) encapsulates the dependence on the scale and is given 
by an integral over the primordial bispectrum. 

To summarise, the cosmological parameter vector in this non- 
Gaussian extension of the ACDM model is 


P - {Gn„Cr8,Wo,Wa,0(,,nb,^0,«S,/NLl • (20) 

We assume /nl = 0 as the fiducial value of the non-Gaussian am¬ 
plitude. 

The level of primordial non-Gaussianity has recently been 
constrained to high precision thanks to Planck CMB data 
dPlanck Collaboration et ^Jl2013h . -4 < /nl :£ 11, for the case of 
a local bispectrum shap^. Bounds from galaxy cluster abundance 
show consistency wi th the Gaussian scenario, -91 :£ /nl 78 
dShandera et al.l2013h . Co nstraints from the d istribution of clusters 
are even less stringent dMana et alj l2013h . The clustering of 
Euclid spectroscopic galaxies alone is expected to restrict the 
allowed non-Gaussian parameter spa c e down to AfNL ~ a few 
I Carbone. Verde & Matures j l20Q8bl : IVerde & Matarresel 1 20091 : 
lFedelietalJl201lh . ^ 


^ The Planck CMB constraints on primordial non-Gaussianity have been 
converted here into the LSS convention. 
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4.3 Parametrise deviation from General Relativity 


We studied another extension to the standard ACDM cosmol¬ 
ogy, based on deviations of the law of gravity from GR. As 
a matter of fact, a number of non-stan dard gravity models 
have been proposed in the literature (e.g., Hu & Sawickil 1 20071 : 
iCapozziello & de Laurenti jl201 ll : lAmendola et alj|2013h in order 
to explain the low-redshift accelerated expansion of the Universe 
without need for the DE fluid. Many of these models give rise to 
modifications of the late-time linear growth of cosmological struc¬ 
ture, which can be parametrised as 


rflnD(fl) 

din a 


= tiLG), 


( 21 ) 


where y is dubbed the growth index (e.g. iLahav et alJfT^ . GR 
predi cts a nearly c onstant and scale-independent value of y 0.55 
(e.g. lLindei|[2005h . Significant deviations from this value would 
hence signal a violation of the standard theory of gravity on large 
scales. The corresponding vector of cosmological parameters in 
this case reads 


{Dm,cr8,Wo,Wa,f2A.,nb,Ho,ns,'Xl , 


( 22 ) 


with y = 0.55 taken as o ur reference value. U sing number counts 
of X-ray clusters alone, iMantz et al.l (1201 5h have found values 
of y consistent with GR (y = 0.48 ± 0.19). From a sample of 
SZ-sel ected clusters in SPT survey y = 0 . 73 ± 0.28 has been 
found dSocQuet et alj|20l3) . iLombriser et al (200) have directly 
constrained the f(R) model bv iHu & Sawicki ( 2007 ) by exploiting 
an optically selected cluster sample. 


4.4 Model with non-minimal neutrino mass 


In our analysis we also consider the case of massive neutrinos, with 
the associated density parameter Qy as the relevant parameter to 
be constrained, fly is related to the total neutrino mass, niyj, 
through the relation: 


Pv 

'pc ~ 93.14h2eV’ 


(23) 


where py and pc are the z = 0 neutrino and critical mass densi¬ 
ties, respectively, and Ny is the number of massive neutrinos. A 
larger value of fly acts on the observed matter power spectrum 
in two ways (e.g. lUesgourgues & Pastorll200 ^ : lMarulli et alJl201ll : 

iMassara, Villaescusa-Navarro & Viej2014h . The peak of the power 

spectrum is shifted to larger scale, because a larger value of the 
radiation density postpones the time of equality. Moreover, since 
neutrinos free-stream over the scale of galaxy clusters, they do not 
contribute to the clustered collapsed mass on such a scale. As a 
consequence, the halo mass function at fixed val ue of fl„, will be 
below the one expected in a purely CDM model. iBrandbyge et alj 
( l2010l) have shown that results from N -body simulations with mas- 
siv e neutrinos can be r eproduced in a more accurate way by using 
the lTinker et alj ( l2008h halo mass function with 


Pm PCDM + Pb - Pm- Pv, (24) 

where Pm,PcDM,Pb and py are the total mass, CDM, and baryon 
and neutrino densities. Based on the analysis of a n extended set 
of N-b ody simulations. ICastorina et al] ( |2014|) and ICostanzi et al] 
( l2013bl) have shown that, since neutrinos play a negligible role in 
the gravitational collapse, only the contribution of cold dark matter 
and baryons to the power spectrum has to be used to compute the 


r.m.s. of the linear matter perturbations, a(R), in the computation 
of the halo mass function and linear bias: 


Pm PcDInik) - P m(k) 


kicDMTcDMik, z) -I- QhTbik,z) 

(Dj -I- Q.CDM)Tm(k, z) 


(25) 


Here Tcdm, Tb and T,„ are the CDM, baryon, and total matter trans¬ 
fer functions respectively, and P„, is the total matter power spec¬ 
trum. 

Hence, the cosmological parameter vector we use in this case 
is: 


p = {Om.o'g.wo, Db.tfoPJs.flvl , (26) 

with a fiducial value of fly = 0.0016 that corresponds to 

y, rUy = 0.06 for t hree degenerate neutrinos dCarbone et al.ll20l^ : 
iMantz et al.l l2015h . Currently, great attention has been devoted 
to derive constraints on the neutrino mass from the combination 
of galaxy clusters with other LSS observables. The analysis of 
the Planck SZ cluster sample resu lt ed in Y^niy = 0.20 ± 0.09 eV 
dPlanck Collaboration et al.ir2014cl) . IMantz et ^ d2014h . combin¬ 
ing cluster, CMB, SNla and BAG data, fo und E Wy < (3 .38 eV 
at 95.4 per cent c.l. in a wCDM universe. ICostanzi et alj d2014ll 
found 2 my < 0.15 eV (68 per cent c.l.) in a ACDM universe, 
for a three active neutrino scenario, usi ng cluster counts, CM B, 
BAG, Lyman-ff, and cosmic shear data. In lBocguet et alj d2014l) the 
analysis of SPT cluster sample resulted in 2 niy = 0.148 ± 0.081 eV. 


4.5 Parameters of the mass-observable scaling relation 

In our FM analysis, besides the cosmological parameter vectors de¬ 
tailed above, we also include four extra parameters to model intrin¬ 
sic scatter and bias in the scaling relation between the observed and 
true galaxy cluster masses (see equation above). We assume the 
following parametrisation for the bias and the scatter, respectively: 


In Mb.as(z) = Sm.o + a In (1 -(- z) 


and 


= + (27) 

We select the following fiducial values 

Pnuisance,F = = 0, tt = 0, = 0.2,^6 = 0.125) . (28) 

We refer to these four parameters as nuisance parameters hence¬ 
forth. With the fiducial nuisance parameter vector there is no bias 
in the true mass-observable r elation and the value of the scatter at 
z = 0 is in accordance with iRvkoff et alj d2012h . Also, the fidu¬ 
cial value for [} makes the scatter increase with redshift, reach¬ 
ing cTbiAf - 0.6 at the maximum redshift of the Euclid survey 

(Zmax = 2). 

Currently, the mass-observable relations are not known over 
the full redshift range that will be covered by Euclid. In the Euclid 
survey it will be possible to calibrate such relation with its uncer¬ 
tainties thanks to the weak lensing and spectroscopic surveys. We 
estimate that Euclid has the potential to calibrate the scaling rela¬ 
tion to < 15 per cent accuracy out to z < 1.5 (see Appendix iBt. 

In the following Section, we will consider the two extreme 
cases where we assume (i) no prior information on the nuisance 
parameters, and (ii) perfect knowledge of the mass-observable re¬ 
lation. 
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Figure 3. Number of clusters above a given redshift to be detected with 
overdensities N field > 5 and > 3 in the Euclid photometric survey 
(dotted blue and solid red lines, respectively). We also show the number 
density of clusters expected to be detected within redshift bins of width 
Aj = 0.1 for the same detection thresholds (dotted cyan and solid magenta 
histograms, respectively).The numbers have been obtained by using the ref¬ 
erence values of cosmological and nuisance parameters (see Section|4). 

5 RESULTS 

Here, we present the constraints on the cosmological parameter 
vectors introduced in the previous Section, using the FM formal¬ 
ism. As a first result, we plot in Fig.l^the histograms corresponding 
to the redshift distributions, n{z) = AzdN/dz (equatiorO, of Euclid 
photometric galaxy clusters, obtained by adopting the two selection 
functions, which correspond to the two different detection thresh¬ 
olds Nioo.clo'fieid > 3 and 5 (see Fig.j^, and by using the reference 
values of cosmological and nuisance parameters. The two curves 
show the corresponding cumulative redshift distributions, n(> z), 
i.e., the total number of clusters detected above a given redshift. 
Euclid will detect ~ 2 x 10^ objects with A^soo.c/crfleid ^ 5 at all 
redshifts, with about ~ 4 x 10'^ of them at z > 1. By lowering the 
detection threshold down to Alsoo.c/crfieid = 3, these numbers rise 
up to ~ 2 X 10® clusters at all redshifts, with ~ 4 x 10® of them 
at z > 1. The large statistics of clusters at z > 1 provides a wide 
redshift leverage over which to follow the growth rate of pertur¬ 
bations. As a comparison, DES will detect ~ 1.7 x 10® clusters 
(with more than 10 bright red-sequence galaxies) and with masses 
greater than ~ 5 x lO'® Mq out to z ~ 1.5 in the survey area of 5000 


deg® 0 eROSITA dPillenich, Porciani & Reipri^l2012h will detect 
~ 9.3 X 10'* clusters with masses greater than ~ 5 x 10*® Mq in the 
survey area of 27.000 deg®, almost all at z < 1. 

In Figs. 016117] and[8]we show the forecasted constraints from 
Euclid photometric clusters on suitable pairs of cosmological pa¬ 
rameters. The ellipses in these figures always correspond to the 
68 per cent c.l. after marginalisation over all other cosmological 
parameters and nuisance parameters. In each of these figures, the 
blue dotted contours are obtained by combining the number counts 
(NC) FM (equation [3]l and the cluster power spectrum (PS) FM 
(equatior®, assuming no prior information on any of the cosmo¬ 
logical and nuisance parameters. Also, the cluster sample is defined 
by the selection Asoo.c/o'flcid S' 3. The green dash-dotted contours 
are obtained in the same way except for the addition of strong priors 
on the nuisance parameters, i.e. assuming perfect knowledge of the 
scaling relation between the true and the observed cluster mass (this 
is labelled as ”-l-known SR” in the figures). The magenta solid con¬ 
tours have been obtained by further introducing prior information 
from Planck data (label-ed ”+Planck prior” in the figures). Finally, 
the cyan solid contours represent the same combination of informa¬ 
tion as the magenta solid ones (NC-t-PS-l-known SR-t- Planck prior) 
obtained from the cluster sample with selection corresponding to 
A^doo.c/ffieid > 5. In the figures, we indicate these contours with the 
label 5cr. 

When using the Planck priors, we take for the CPL DE model 
the correlation matri x obtained by combining P l anck C MB data 
with the BAOs from IPlanck Collaboration et aP ( l2014aH^ for the 
parameters of the ACDM cosmology (assuming fl* = 0), plus 
Wo and vvj^- For the non-Gaussian case, we use priors from the 
Planck obtained for the ACDM model plus parameter^ We 
also added a flat prior on the level of non-Gaussianity correspond¬ 
ing to -5.8 < ^ Finally^ for the modified gravity and 

the neutrino scenario we also used priors from the Planck analysis 
carried out over the parameters of the ACDM model plus 

In Fig.0 we show the constraints on Dm and erg (left panel), 
as well as those on the two CPL DE parameters wq and Wa (right 
panel). The contours on the - <Tg plane for the combination of 
number counts and clustering of Atsoo.c/cTfieid > 3 galaxy clusters 
are rather tight. The information provided by the number density 
of clusters alone defines the degeneracy direction between and 
(Tg, with the following constraints: AD,„ = 0.009, Acrg = 0.006. 
Information from the cluster power spectrum alone does not pro¬ 
vide stringent constraints on the - erg plane. However, using the 
combination of the PS with NC FM, the values of both parameters 
are constrained to high accuracy: Afl„, = 0.0019, Acrg = 0.0032 
(see Table[T). By assuming a perfect knowledge of the scaling re¬ 
lation between true and observed cluster mass, the bounds improve 
significantly. This is especially true for erg, which is more affected 
by the nuisance parameters than Dm- Including information from 
the Planck priors does not improve the forecasted constraints sig¬ 
nificantly, in keeping with the expectation that the Euclid cluster 
bounds are, by themselves, competitive with CMB bounds. 

Taking the ACDM as a reference model, its parameters will 
be constrained with a precision of ~ 10“®, 

AD,„ = 5.9 10^^ Acrg = 4.9 10-^ Ah = 7.2 10^\ 


® https://www.darkenergysurvey.org/reports/proposal-standalone.pdf 

® Available at http://wiki.cosmos. esa. int/planckpla/index.php/Cosmological_P 
® PLA/base_w_wa/planck_lowLlowLike_BAO 
* PLA/base_omegak/planck_lowlJowLike 
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Figure 4. Constraints at the 68 per cent c.l. on the parameters Qm and erg (left panel) and on the parameters wq and Wa for the DE EoS evolution (right 
panel). In each panel, we show forecasts for the A^ 5 oo,c/(Afield ^ 3 Euclid photometric cluster selection obtained by (i) NC, the EM number counts (red 
dash-dotted contours), (ii) NC+PS, the combination of EM NC and power spectrum (PS) information (blue dotted contours), (hi) NC+PS+known SR, i.e. by 
additionally assuming a perfect knowledge of the nuisance parameters (green dash-dotted contours), and (iv) NC+PS+known SR+Planck prior, i.e. by also 
adding information from Planck CMB data (magenta solid contours). With cyan solid lines we show forecasts for the A^ 5 oo,c/<^fieid ^ 5 Euclid photometric 
cluster selection in the case NC+PS+known SP.+Planck prior (labelled 5 (t). Planck information includes prior on ACDM parameters and the DE EoS 
parameters. 


An* = 8.4 10^^ An, = 3.3 10“^ (29) 

thanks to the unprecedented number of clusters that will be de¬ 
tected at high redshift. These constraints have been obtained with 
the A^soc.r/cfieid > 3 selection function, from cluster number counts 
and power spectrum, by assuming strong prior on the nuisance pa¬ 
rameters, and no prior from Planck. 

These results emphasise the importance of exploring the high- 
redshift clusters in survey mode. Of course a good knowledge of the 
astrophysical process taking place in clusters is fundamental to cali¬ 
brate the mass-observable scaling relations, and also to optimise the 
detection algorithms. Hence detailed follow-ups of restricted sam- 


pies 

of clusters (such as, e.g., CLASH, CCCP, WtGiPostman et al, 

2012 

;lRosati et alj2014l;lHoekstra et alj2012l;lvon der Linden et al 

2014 

) retain a crucial importance. 


On the other hand, the inclusion of Planck priors shall bring 
a substantial improvement over the bounds to the DE parameters. 
This result is expected, since the CMB data provides stringent con¬ 
straints on the curvature, thereby breaki ng the degeneracy be tween 
n* and the evolution of the DE EoS jSartoris et al.ll201^ . The 
contribution of the PS information is less important for (wo,Wa) 
with respect to (Dni.crs): however, the FoM increases from ~ 30 
in case of NC alone to ~ 73 for NC-hPS constraints (see Table[T]l. 
For both DE EoS parameters, it is crucial to have a well calibrated 
scali ng relation over the redshift range sampled by the cluster sur¬ 
vey dSartoris et al.ll201^ . Indeed, by combining NC and PS, and 
assuming perfect knowledge of the scaling relation increases the 
FoM to - 291. When we include the Planck data, i.e. we set a prior 
on the curvature, we obtain FoM= 802, with Awo = 0.017 and 
AWa - 0.074 (see Table[TJ. 

When we restrict our analysis to the wCDM model (that is 


characterised by the six free parameters fl,„, trg, A, fl*, n,, w), we 
obtain Aw = 0.005. If we also add as a free parameter, we ob¬ 
tain Awo = 0.013 and Aw^ = 0.048. These constraints have been 
obtained with the Vsoo.c/o'fleid > 3 selection function, from cluster 
number counts and power spectrum, by assuming strong prior on 
the nuisance parameters, and no prior from Planck. 

In both panels of Fig. [4] the adoption of a more conservative 
cluster selection (A^soo.c/crfleid > 5) significantly worsens the fore¬ 
casted cosmological constraints. For instance, the FoM is degraded 
down to 209 in the best-case scenario, as a consequence of the sig¬ 
nificantly degraded statistics corresponding to the higher selection 
threshold. 

In Fig.|5] we show how the FoM depends on the limiting red¬ 
shift of the survey. The FoM shown in this figure refers to number 
counts (NC) in the Njoo.c/crfleid > 3 Euclid photometric cluster se¬ 
lection. The FoM for a survey reaching out to z < 1.2 is only half 
the FoM of an equivalent survey reaching out to z < 2. It is there¬ 
fore important that the redshift range covered by the survey be large 
enough to allow a comparison of the behaviour of DE over a suffi¬ 
ciently long cosmological timescale. In this sense, the Euclid sur¬ 
vey will have a unique advantage over other existing and planned 
surveys. 

In Fig. we show cosmological constraints in the expanded 
parameter space which includes non-Gaussian primordial density 
fluctuations. Specifically, we display the constraints in the /ml - Cg 
plane. Thanks to the peculiar scale-dependence that primordial 
non-Gaussianity induces on the linear bias parameter, the power 
spectrum of the cluster distribution turns out to be much more sen¬ 
sitive to /ml than it is to trg. This is clearly demonstrated by the red 
dash-dotted contour, which shows forecasted constraints derived 
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Table 1. Figure of Merit (FoM) and constraints on cosmological parameters as obtained by progressively adding the FM information for different models, 
for two different detection thresholds (A^ 5 oo,c/<^fieid ^ 3 and 5). Constraints are shown at 68 per cent c.l. after marginalisation over all other cosmological 
parameters and nuisance parameters in the aiTays. 


A^ 50 o,c/<^fieid ^ 3 Euclid photometric cluster selection 

Parameter an'ays: 


Eqs.[T6]&[28] 

Eqs.[22]&|28] 

Eqs.[20]&|28] 

Eqs.[26]&|28] 

Constraints: 

FoM 

Awo AWfl 


Ao-g 

Ay 

A/jvl 

AHv 

NC+PS 

73 

0.037 0.38 

0.0019 

0.0032 

0.023 

6.67 

0.0015 

NC+PS+known SR 

291 

0.034 0.16 

0.0011 

0.0014 

0.020 

6.58 

0.0013 

NC+PS+known SR+Planck 

802 

0.017 0.074 

0.0010 

0.0012 

0.015 

4.93 

0.0012 

A^ 50 o,c/<^fieid ^ 5 Euclid photometric cluster selection 

NC+PS+known SR+Planck 

209 

0.034 0.12 

0.0022 

0.0026 

0.034 

6.74 

0.0020 


1 I-1-1-1-1-1-1-1- - 

0.9 - • - 

0.8 - • - 

0.7 - 

rj • 

II 

I 0.6 - • 

N 

S 

o 0.5 - * 

H- 

2 0.4 - 

o 

n- 

0.3 - 

0.2 - * 

0.1 - 

0 _ S _I_I_I_I_I_I_I_ 

0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

^max 


Figure 5. Relative FoM for number counts in the Wjoo.c/o'fleld ^ 3 Euclid 
photometric cluster selection, as a function of the limiting redshift z^ax of 
the survey, i.e. the ratio between the FoM evaluated over 0.2 < z < Zmai and 
the FoM evaluated over 0.2 < z < 2.0. 

from cluster clustering alone. Quite clearly, erg is basically uncon¬ 
strained on the scale of the figure, while /nl is constrained with 
an uncertainty A/nl ~ 7.4. The addition of cluster number counts 
changes very little the bounds for primordial non-Gaussianity, how¬ 
ever it improves substantially those for the amplitude of the mat¬ 
ter power spectrum (see Table [TJ. This helps to define the degen¬ 
eracy between /nl and erg that are both related to the timing of 
structure formation. Interestingly, the estimation of primordial non- 
Gaussianity is weakly sensitive to the nuisance parameters. Indeed, 
when a perfect knowledge of the scaling relation between true and 
observed cluster mass is assumed, only the constraints on trg im¬ 
prove significantly. Planck priors does not affect substantially the 
constraints on /ml. 

When we restrict our analysis to the ACDM model plus the 
non-Gaussianity parameter /nl, we obtain A/nl = 6.44. This con¬ 
straint has been obtained with the V 5 oo,c/cr(i;.Li > 3 selection func¬ 
tion, from cluster number counts and power spectrum, by assuming 
strong prior on the nuisance parameters, and no prior from Planck. 



Figure 6. Constraints at the 68 per cent c.l. on the /nl - erg parameters. We 
show forecasts for the Asoo.c/crfieij > 3 Euclid photometric cluster selec¬ 
tion obtained by (i) PS, the FM power spectrum (red dash-dotted contours), 
(ii) NC+PS, the combination of FM number counts (NC) and PS informa¬ 
tion (blue dotted contours), (iii) NC+PS+known SR, i.e. by additionally as¬ 
suming a perfect knowledge of the nuisance parameters (green dash-dotted 
contours), and (iv) NC+PS+known SR+Planck prior, i.e. by also adding 
information from Planck CMB data (magenta solid contours). With cyan 
solid lines we show forecasts for the Wjoo.c/<rfield > 5 Euclid photomet¬ 
ric cluster selection in the case NC+PS+known SR+P!anck prior (labelled 
5cr). Planck information includes prior on ACDM+f2t+/NL parameters. 

Forecast for eROSITA dPillepich, Porciani & ReiDri^l2012h pre¬ 
dict a similar precision, since the narrower redshift range of this 
survey (with respect to Euclid) is compensated by its wider area, 
which allows a better sampling of large scale modes. 

We point out that in this analysis we are assuming the most 
commonly used parametrisation of non-Gaussianity, where /nl is 
considered scale-invariant. However, there are models that predict 
otherwise. For these, the combination of clusters and CMB data 
complement each other well, providing tight constraints on the pos¬ 
sible scale dependence of /nl- 

As for the models including GR violation, we show in Fig. 
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Figure 7. Constraints at the 68 per cent c.l. on the y - erg parameter 
plane. We show forecasts for the Alsoo.c/cfieid > 3 Euclid photometric 
cluster selection obtained by (i) NC+PS, the combination of FM number 
counts (NC) and power spectrum (PS) information (blue dotted contours), 
(ii) NC+PS+known SR, i.e. by additionally assuming a perfect knowl¬ 
edge of the nuisance parameters (green dash-dotted contours), and (iii) 
NC-rPS-rknown SR+Planck prior, i.e. by also adding information from 
Planck CMB data (magenta solid contours). With cyan solid lines we show 
forecasts for the A^soo.c/o'field > 5 Euclid photometric cluster selection in 
the case NC-l-PS-l-known SR+Planck prior (labelled 5cr). Planck informa¬ 
tion includes prior on ACDM-rfl^t parameters. 


|7]the constraints on erg and the growth parameter y. Similarly to 
the constraints on the n^-erg plane, the constraints on y are not 
strongly affected by the inclusion of Planck priors, thus implying 
that galaxy clusters are by themselves excellent tools to detect sig¬ 
nature of modified gravity through its effect on the growth of pertur¬ 
bations. Significant degradation of the constraining power happens 
if a higher threshold for cluster detection is chosen. 

Restricting our analysis to the ACDM model plus the y pa¬ 
rameter we obtain Ay = 0.006. This constraint has been obtained 
with the Asoo.c/cTfldd > 3 selection function, from cluster number 
counts and power spectrum, by assuming strong prior on the nui¬ 
sance parameters, and no prior from Planck. 

Finally, we show in Fig. the joint cosmological constraints 
on (Tg and the neutrino density parameter fly. The presence of neu¬ 
trinos with masses in the sub-eV range requires higher values of 
(Tg: increasing at fixed has the effect of shifting the epoch 
of matter-radiation equality to a later time and to reduce the growth 
of density perturbations at small scales in the post-recombination 
epoch. As a consequence, a larger value of erg is required to com¬ 
pensate these effects. We use the Planck prior mainly to add infor¬ 
mation on the geometry of the Universe, and the standard ACDM 
parameters. We obtain the constraints AQy = 0.0012 (correspond¬ 
ing to AJjnjv - 0.01 ). The constraints on the neutrino density 
parameter are weakly affected by the inclusion of a prior on the 
nuisance parameters. Flowever, there is a degradation by a factor 
of ~ 2 of the constraining power if the selection function with the 
higher threshold for cluster detection is chosen (see Table[T}. 





Figure 8. Constraints at the 68 per cent c.l. in the fly — erg parameter 
plane. We show forecasts for the A^soo.c/Tfield > 3 Euclid photometric 
cluster selection obtained by (i) NC-l-PS, the combination of FM number 
counts (NC) and power spectrum (PS) information (blue dotted contours), 
(ii) NC-l-PS-l-known SR, i.e. by additionally assuming a perfect knowl¬ 
edge of the nuisance parameters (green dash-dotted contours), and (iii) 
NC-l-PS-l-known SR+Planck prior, i.e. by also adding information from 
Planck CMB data (magenta solid contours). With cyan solid lines we show 
forecasts for the A^soo.c/ffieid > 5 Euclid photometric cluster selection in 
the case NC-l-PS-l-known SR+Planck prior (labelled 5cr). Planck informa¬ 
tion includes prior on ACDMd-Dt parameters. 

To gauge the impact of a particular choice of the selection 
function on the cosmological constraints, we have so far shown our 
results for both the Ajoo.c/crflcid > 3 and the Ajoo.c/crfleid > 5 Eu¬ 
clid photometric cluster selection functions. As a further test, we 
consider the etfect on the (h'o,Wo) constraints of adopting a flat 
selection function with log(M 2 oo,c) = 13.9, within 0.2 < z < 2. 
With this flat selection function there are less clusters than with the 
^soo.c/o'fieid > 3 one, both in total (~ 1.4 x 10* vs. ~ 1.6 x 10*, 
respectively) and within 0.4 < z < 1.2. However, the number of 
clusters at z > 1 is higher (~ 3.2 x 10*) for the^af selection func¬ 
tion than for the Asoo.c/o'fieid > 3 one (~ 1.9 x 10*). The effect of a 
larger number of high-z clusters in the flat selection function sam¬ 
ple compensates for the smaller total number of clusters in provid¬ 
ing similar constraints on the cosmological parameters to those ob¬ 
tained with the Asoo.c/cTfieid > 3 selection function sample (changes 
are < 10% on the constraints on the DE parameters). This suggests 
that the precise shape of the selection function has little impact on 
our results, and in any case much less than its overall normalisation. 


6 CONCLUSIONS 

In this paper, we presented a comprehensive analysis of the fore¬ 
casts on the parameters that describe different extensions of the 
standard ACDM model. These were based on the selection func¬ 
tion of galaxy clusters from the wide photometric survey to be car¬ 
ried out with the Euclid satellite, a medium-size ESA mission to 
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be launched in 2020. We presented the derivation of this selection 
function and the Fisher Matrix formalism employed to derive cos¬ 
mological constraints. This is the same formalism that has been 
used in the Euclid Red Book iLaureiis etalj|201lh . Our main re¬ 
sults can be summarised as follows. 


• Using photometric selection, we found that Euclid will detect 
galaxy clusters at A^soo.c/crfleid > 3 with a minimum mass of ~ 0.9 - 
1 X IO'^^Mq. As a result, the Euclid photometric cluster catalogue 
should include ~ 2 x 10^ objects, with about one fifth of them at 

z ^ 1. 

• The Euclid cluster catalogue has the potential of providing 
tight constraints on a number of cosmological parameters, such as 
the normalisation of the matter power spectrum erg, the total matter 
density parameter flni, a redshift-dependent DE equation of state, 
primordial non-Gaussianity, modified gravity, and neutrino masses 
(see Table[Tll. We predict that most of these constraints will he even 
tighter than current bounds available from Planck. The constrain¬ 
ing power of the Euclid cluster catalogue relies on its unique broad 
redshift coverage, reaching out to z = 2. 

• Knowledge of the scaling relation between the true and the ob¬ 
served cluster mass turns out to be one of the most important factors 
determining the constraining power of the Euclid cluster catalogue 
for cosmology. The Euclid mission will have a distinct advantage 
in this respect, namely the possibility to calibrate such relation, at 
least up to z = 1.5, with < 10 and < 30 per cent accuracy, us¬ 
ing the weak lensing and spectroscopic surveys, respectively (see 
Appendix iBt. The deep Euclid survey will allow to extend the cali¬ 
bration to even higher redshifts, although with lower precision than 
in the wide survey, due to lower number statistics. 


With the future large surveys, like Euclid, that will be car¬ 
ried out with the next generation of telescopes, the number of de¬ 
tected clusters from the individual surveys will range from thou¬ 
sands to tens of thousands. As we have shown in this paper, this 
will allow to constrain most cosmological parameters to a preci¬ 
sion level of a few per cent. Currently, theoretical halo mass func¬ 
tions are defined wi th an uncertainty o f ~ 5 per cent in the standard 
ACDM model (e.g. [Tinker et al.l2008l) . and many efforts have been 
devoted in the last years to better sample the high mass regime 
dWatson et al.ll201^ . To maximally extract cosmological informa¬ 
tion from these cluster surveys, it becomes critical to specify the 
theoretical halo mass function to better than a few percent accu¬ 
racy for a range of cosmol ogies. A substantia l effort is currently 
ongoing in this direction l Grossi_etaO 120071; [DaMetaQ |2()oj; 


Cui jB^di & Borganill2012l : Lombriser et alj|2013l : Castorina et alj 


2014h . Moreover, cosmological hydrodynamic simulations will 
have to precise the impact of baryons on the shape of the mass 
profile, which has already been shown to be quite significant 
(iRudd. Zentner & Kravtso^ l2008l: IStanek. Rudd & Evrard 1 20091 : 

ICui et alj|2012l: ^i. Borgani & Murantell2014). 
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Figure Al. Number of cluster galaxies with spectroscopic redshifts within 
i'200,c expected in the Euclid survey, as a function of redshift for clusters 
of different masses, log(M 2 oo,c/JJQ) = 15.0,14.5,14.0,13.5, solid, dot- 
dashed, dashed, dotted lines, respectively. These numbers are for the case 
of an evolving Ho- luminosity function also beyond z= 1.3, i.e. they corre¬ 
spond to the solid blue curve in Fig. lA2l ftoD panel). 



APPENDIX A: THE EUCLID SPECTROSCOPIC SURVEY 


We use a procedure similar to the one described in Sectionj^to de¬ 
termine the number of spectroscopic cluster galaxies within r 2 oo.c, 
as a function of both cluster mass and redshift. Since the Euclid 
spectroscopic survey is flux-limited in the Her line, we consider the 
cluster Ha LF. There are not many determ inations of the cluster 
Ha LF in the literature. We use the resultsj)f[l glesms-Paramo et alj 
(I 2 OO 2 L for two nearby clusters, z = 0.02), Balogh et alJ 1 20021 . for a 
z = 0.18 rich clust er), Umeda et"^ (I 2 OO 4 I for a z = 0.25 cluster), 
and iKodama et al.l (|2004 for a z = 0.4 cluster). 

The redshift evolution of the cluster Ha LF is (at best) poorly 
constrained, hence we have to make several assumptions for its 
three parameters, the characteristic luminosity L*, the normalisa¬ 
tion 0*, and the faint-end slope a. We consider two possible evolu¬ 
tions. In the first, we assume the z-evolution of L* to be the same 


as the one measured for the field Ha LF, i.e. L* oc ( 1 + for 
z < 1.3, and no further evolution at higher redshift jGeach et alj 
I 2 OICI) . In the second, we allow L* to evolve at z > 1.3 with the 
same z-dependence established at lower redshifts. The second sce¬ 
nario is based on the idea that the preferred sites for galaxy star- 


Figure A2. Selection function for the Euclid spectroscopic survey. In the 
top panel the solid blue curve indicates the selection function for clusters 
with 5 galaxies with measured spectroscopic redshift within r 2 oo,c- This 
curve depends on the assumption that L* continues to evolve bey o nd z = 1 -3 
following the same evolution law determined by iGeach et alJ bOlOl) for 
lower z. The dash-dotted curve depends instead on the assumption that there 
is no further evolution of L * beyond z = 1-3, consistently with what is ob¬ 
served for the field Ho- LF iGeach et aflllOlOl) . The dashed red curve is an 
independent estimate based on the the number densities of Ho' field galaxies 
per redshift bin, estimated by Pozzetti et al. (in prep.). In the lower panel 
the solid, dash-dotted and dashed lines show results for clusters with at least 
5, 10, and 20 galaxies, respectively, with measured spectroscopic redshift 
within r 2 oo,c, based on the assumption that L* continues to evolve beyond 
z = 1.3. The dotted line is the selection function for the Euclid photometric 
survey (from Fig. [3, shown as a reference. 
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formation tends to shift to higher-density regions at higher redshifts 
telbaz et alj|2007ti . e ven if the redshift a t which this shift occurs is 
not well constrained IZinaro et alll2014ll 

The different cluster LFs we consider have been determined 
for different overdensities, A. To evaluate the A = 200 value of L* 
at z = 0, we perform a regression analysis between log [L*/(l -l- 
Zc)^ *] and log A. We find = 3.8 x lO'** erg s“'. Similarly to 
what we did in Section for the Ks LF, we assume (p* oc tC{z). 
We then take the average of the 0* values obtained for the different 
clusters, after rescaling them for the factor 200//o/[A//(z)], and 
find = 1.1 Mpc“^. As for a, we fix it to the value -0.7 ob tained 
for the two nearby clusters bv Ilglesias-Paramo et all l l2002h . since 
the other clusters observations were not deep enough to constrain 
the FItt LF faint-end. 

We convert the Tier luminosities into fluxes using /hq. = 
f'Ha/(2 X AnD^yj, where is the cluster luminosity distance 
and the factor 1/2 accounts for the average dust extinction 
jKodama et alj|2003l . By integrating the LF down to the flux limit 
of the Euclid spectroscopic survey (3 x 10“'® erg s“' cm“-), we fi¬ 
nally obtain the expected number density of galaxies within r2oo,c- 
By multiplying the number density of galaxies within rioo.c by the 
volume of the sphere of radius r2oo,c, we obtain the number of galax¬ 
ies in a cluster with Flor flux above the Euclid survey limit. Finally, 
we multiply this number by the expected completeness of the spec¬ 
troscopic survey, 80 per cent. 

In Fig. lAll we show the resulting estimates of the number of 
cluster galaxies with spectroscopic redshifts within r2oo,c, as a func¬ 
tion of redshift for clusters of different masses, for the case of an 
evolving Her LF beyond z = 1.3. Note that only the redshift range 
0.9-1.8 is shown, since this is the detectability range of the Hcr-line 
in the Euclid survey, according to the current design baselinfl 

We also consider the following, independent estimate of the 
cluster selection function in the Euclid spectroscopic survey. We 
use Pozzetti et al.’s (in prep.) estimates of the number densities 
of Ho-emitting field galaxies per square degree and redshift bin, 
that we convert to volume densities, Ufj. To estimate the expected 
number density of Hcr-emitting galaxies in a cluster, we used 
«£■/ = nfdb(z)Apclp,„, where pc is the critical density and p„, the 
mass density of the Universe at any given redshift, A is the over¬ 
density we want to sample in the cluster, and b(z) is the redshift- 
dependent bias parameter that accounts for the different distribu¬ 
tion of Her galaxies and the underlying matter distribution. Taking 
A = 200, the number of Her galaxies in a cluster of mass M 2 Q 0 .C 
is A(< r2oo,c) = (AzT/3)rlgQ^nci. We estimate the bias b(z) from the 
comparison of the real-space correlation functions of matter and 
Her galaxies, b = where 7 is the slope of the correla¬ 

tion function. We use the correlation lengths of the diffuse matter 
in our adopted cosmology, and those of Her galaxies with lumi¬ 
nosities cor responding to the E uclid flux limit at any given redshift 
(taken from ISobral et aklboioh . We estimate b{z = 0.9) = 1.9 and 
b{z = 2.0) = 3.5, and interpolate b{z) between these two values at 
any redshift in the range 0.9-2.0. 

In Fig. lA2l we show the limiting mass M200.C of a cluster with 
at least galaxies with measured spectroscopic redshift within 
''200,c as a function of the cluster redshift. This is the selection func¬ 
tion of clusters in the Euclid spectroscopic survey, in the sense 
that N, concordant redshifts within a region of typical cluster size 


^ See the "Euclid GC Interim Science Review” by Guzzo & Percival, at 
http://internal.euclid-ec.org/?page_id=714 Access restricted to the Euclid 
Consortium members. 



Figure Bl. Constraints at the 68 per cent c.l. in the Wo — wq parameter 
plane. We show forecasts for the Asoo.c/crjieid > 3 Euclid photometric clus¬ 
ter survey obtained by (i) combining the FM information for number counts 
and power spectrum (NC+PS; blue dotted contour), (ii) same as (i) but as¬ 
suming perfect knowledge of the evolution of the scatter (see equation 1271 
orange dashed contour); (hi) same as (i) but assuming perfect knowledge 
of the evolution of both the scatter and the bias (black solid contour); (iv) 
same as (i) but assuming perfect knowledge of all the four nuisance param¬ 
eters (green dash-dotted contour). The blue and green curves are the same 
of Fig.|4]right panel. Note that the solid black and the green dashed ellipses 
are almost coincident. 

(i.e., r 2 oo,c) are required to identify a cluster. The three different es¬ 
timates of the spectroscopic selection function for clusters in the 
Euclid survey are rather different, and this reflects the current sys¬ 
tematic uncertainties. From Fig. IA2I (bottom panel), one can see 
that the spectroscopic survey selection function is above the photo¬ 
metric survey selection function. Hence, it will prove less efficient 
to search for clusters in the Euclid spectroscopic survey than in the 
photometric survey. Data from the spectroscopic survey will still be 
useful to confirm clusters detected in the photometric survey, thus 
improving the reliability of the sample. 


APPENDIX B: CLUSTER MASS CALIBRATION 

The impact of nuisance parameters on cosmological constraints 
from Euclid photometric clusters is going to be quite significant. 
This is especially true for the parameters directly related to the 
growth of structure history like the matter power spectrum normal¬ 
isation CTg, and for the CLP DE parameter w^, that is particularly 
sensitive to the level of knowledge of the scaling relation evolu¬ 
tion. In Fig. im we show how the cosmological constraints on the 
DE equation of state depend on our knowledge of the scaling re¬ 
lation. In particular, we show that strong constraints on the evolu¬ 
tion of the scatter and the mass bias, allow to greatly improve the 
constraints on the DE EoS parameters. On the other hand, precise 
knowledge of these parameters at z = 0 is not of crucial impor¬ 
tance, as shown by the overlapping constraints in the Wo,Wa plane 
in the figure (solid black and dashed green ellipses). 
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To maximise the scientific return of the Euclid galaxy clus¬ 
ter catalogue, it is therefore very important to know the mass scal¬ 
ing relation in an as much as possible precise and unbiased way. 
There are two avenues to obtain this goal. The first one is to cross¬ 
correlate the Euclid cluster sample with samples obtained at dif¬ 
ferent wavelengths by different surveys. For instance, by the time 
Euclid will fly, the eRosita full-sky X-ray cluster catalogue will be 
available, and will provide an important contribution to the cluster 
true mass estimation. Other useful cluster catalogues will include 
the SZ samples provided by the South-Pole Telescope (SPT), the 
Atacama Cosmology Telescope (ACT), and Planck. 

The second avenue, that represents the strength of the Euclid 
mission, consists in exploiting internal Euclid data. Many photo¬ 
metrically selected clusters will appear as signal-to-noise peaks in 
the Euclid full-sky cosmic shear maps. This weak gravitational 
lensing signal will permit us to estimate the cluster masses with¬ 
out relying on assumptions about dynamical equilibrium. Although 
only the more massive systems will permit individual mass mea¬ 
surements, we can nevertheless statistically calibrate the normali¬ 
sation of the cluster scaling relations down to the lowest masses in 
the catalogue by stacking. An example is given in Fig. IB2I show¬ 
ing the level of precision expected on the mean mass of stacked 
clusters. 

We first measure the mass of individual clusters with a 
matched filter, assuming that the mass density profile of all clus¬ 
ters follows an NFW profile. We then calculate the uncertainty on 
the mean mass of the individual measurements in bins of mass 
(AlogAf 2 oo,c = 0.2) and redshift (Az = 0.1). This result depends 
on the number of clusters expected in e ach bin, and for this purpose 
we hav e adopted the Planck cosmology jPlanck Collaboration et alj 
l2014ah and a Euclid survey of 15,000 square degrees. The figure 
only accounts for shape-noise, with a = 0.3. 

The three curves trace the precision on the mean for mass bins 
centred at Mzoo.^ = 3 x 10 ‘''Mq, 2 x lO'^'Mo, and 1.5 x lO'^'Me 
(from top to bottom) as a function of redshift. We do better on the 
lower mass systems because their larger number compensates for 
their lower individual signal-to-noise measurements. The figure 
demonstrates that Euclid has the potential to calibrate the mean 
mass, and hence scaling relations, to 1% out to redshift unity, and 
to 10% out to z < 1.6 for clusters of M 2 ao,c = 1-5 x lO'"^ Mq. 

At the same time, the spectroscopic part of the Euclid survey 
will provide velocities for a few cluster members in each cluster 
detected with photometric data. Stacking these velocities for many 
clusters in bins of richness and redshift will allow a precise calibra¬ 
tion of the velocity dispersion vs. richness relation, and from this 
of the mass-richness relation. 

In Fig. lB3l we show the number of spectroscopic cluster mem¬ 
bers that will be available for stacks of clusters of given mass in bins 
of Az = 0.1 and A log Mioof — 0.2 (even if, in reality, the stacking 
procedure will be based on mass proxies, such as richness). These 
numbers are evaluated using the spectroscopic selection function 
(bottom panel of Fig. IA2I . and the expected number of clusters 
above a given mass in our adopted cosmology, by considering only 
clusters with at least 5 members with redshifts. In the figure we 
show the predictions for three cluster masses, logM 2 oo,c/^o = 
14.2,14.4,14.6. The curve for logM 2 oo,c/Mo = 14.2 is limited to 
z < 1.25 because of our choice of considering only clusters with 
N, > 5. Note that the curve for log M2oo,c/44o = 14.0 (not shown) 
would be limited to z < 1 (and it would be lie in between those for 
14.2 and 14.4). 

From the analysis of iBiviano et ^ ( l2006h we find that the 
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Figure B2. Calibrating cluster masses with gravitational shear. The curves 
show the expected precision on the mean mass of clusters in bins of 
AlogM 2 oo,c = 0.2 and Az = 0.1, centred on masses (from top to bot¬ 
tom) of M 2 oo,c = 3 x 10*^ Mq (green curve), 2 x IO'^^Mq (red), and 
1.5 x IO'^^Mq (blue). We assume a lensing suiwey of 15,000 sq. deg.^, 
the Tinker mass function in the base ACDM Planck-cosmology, and shape 
noise with cr = 0.3 

statistical noise in the velocity dispersion estimate of a sample 
of ~ 500 cluster members is ~ 9 per cent, which translates into 
a ~ 27 per cent statistical noise in the mass estima t e. A s imilar 
figure has been obtained bv iMamon. Biviano & Boui ll201.3h when 
using the full velocity distribution to constrain cluster masses. 
The value of 500 is displayed in Fig. IB3I and it shows that a very 
precise spectroscopic calibration of cluster masses will be possible 
for stacks of clusters with 14.2 < log M2 oo,c/34g ^ 14.6 over the 
redshift range 0.9 < z 1.2, and even beyond that (z < 1.5) 
for clusters with masses logM 2 oo,clMa - 14.4. Spectroscopic 
calibration of cluster masses at higher redshifts will be feasible 
with reduced precision, but lack of statistics will hamper cluster 
mass calibration at log M2oo,c/4fo < 14.2. 

The wide Euclid survey will allow precise calibration of the 
mass-observable relation out to z ^ 1.6, using gravitational lensing 
and spectroscopy. The deep Euclid survey will allow to extend this 
calibration to even higher redshifts, although with a much more 
limited statistics on the number of clusters. Overall, by combin¬ 
ing Euclid internal mass calibration with the cross correlation with 
external SZ and X-ray surveys, we should be able to significantly 
mitigate the degrading effect of the nuisance parameters on cosmo¬ 
logical constraints. 
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Figure B3. Calibrating cluster masses with spectroscopy. The curves show 
the number of cluster galaxies with redshifts available in stacks of clusters 
in bins of A log M = 0.2 and Az = 0.1, as a function of redshift, for central 
values of the mass bins of log M 2 oo,cIMq = 14.2,14.4,14.6 (red, blue, green 
curves, respectively). The estimate is done only for clusters with a mass 
limit above that required for a minimum of 5 members with redshift - see 
Fig. lA2l bottom panel. This requirement restricts the curve for logM/M© = 
14.2 to z < 1.25. The dotted line shows the value of 500 galaxies as a 
reference. 
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